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The vacuum energy of a conformally coupled scalar field on the d— 
dimensional sphere is calculated. On even spheres it is zero and on 
odd spheres it oscillates in sign. Results for the <i-torus and <i-cube 
are also given. 



1. INTRODUCTION. 



The vacuum energies in ultrastatic space-times T x M have been evaluated for 
various spatial sections and for various fields. A summary would serve no purpose 
here. 

In this paper we wish to present the results of calculations performed some time 
ago but which might be of interest now. Specifically M will be the <i-dimensional 
sphere S d . Results are also given for the <i-cube and <i-torus. 

As usual the field considered is a conformally coupled scalar one satisfying the 
equation, 

with 

e = (d-i)/4d (i) 

and where A 2 is the Laplace-Beltrami operator on M. 

It is shown in Dowker and Kennedy (1978), for example, that the total vacuum 
energy E d is given in terms of the traced zeta function Cd( s ) °f the spatial operator 
-(A 2 - by 

^ = 5& L-*- 1 ^- §)■ P) 

If 1/2) is finite there are no divergences and we shall simply take ^( d (—l/2) 
as the 'renormalised' vacuum energy. Evaluation of the zeta function is now called 
for. 



2. ZETA FUNCTION ON SPHERES. 

This should be a standard problem but the results do not seem to be well 
known. For completeness they are presented here. 

For the choice (1) of £ the eigenvalues A n of — (A 2 — £i?) are 

A n = (n + w- ^) 2 a" 2 , n = 0,l,... (3) 

where 2u = d and a is the radius of S d . The degeneracies are 

r(2w + n-l) /rt 

d ^WTWW) [2ul + 2n - l) 



so that the traced zeta function is, for d > 1, 
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It is necessary to treat odd and even dimensional spheres separately. For odd 
spheres (2ui — l)/2 = q is an integer and the summation can be shifted by setting 
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The second summation clearly gives zero. 

The numerator of the summand is an odd polynomial in n', 



[n' 2 -(q- l) 2 ] . . . [n' 2 - l]n' = J]n ,2r - 1 (-l) r+ X(?) 



which defines the coefficients (c/Cahn and Wolf 1976). 

Thus Cd( s ) can be written as a sum of Riemann zeta functions, 

2a 2s q 



(4) 



Cr(^) has a pole at z/ = 1 of unit residue and no other singularity. Thus Cd( s ) has 
poles at s = r + 1/2, r = 1, 2, . . . , q, with residues 

,2r+l 
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and this is all. 

However according to the general results of Minakshisundaram and Pleijel 
(1949), Cd( s ) should have poles at s = q + 1/2 — n, (n = 0, 1, . . .) with residues 
proportional to a n , the famous expansion coefficients of the heat kernel. What has 
been shown therefore is that for odd spheres with conformal coupling all the a n are 
zero after n = q — 1. This can be confirmed by a direct calculation (G.Kennedy, 
Ph.D. Thesis, University of Manchester, 1979). 

For even spheres on setting d = 2q + 2 and m = 2n + 2q + 1), Cd( s ) can be 
rearranged 
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Again the second sum is zero and the numerator is an even polynomial in m 
of degree 2q, 

[m 2 - (2q - l) 2 ] [m 2 - (2q - 3) 2 ] . . . [m 2 - 1] = £ {-l)* +r 4?- r tf(q) . 

r=0 

This yields 

Cd(s) = (2^! £ C « (2S " 2r " 1 )(" 1 ) r+94_r C 1 " 2 2r " 2S+1 )^) (5) 

as the conformal zeta function on even spheres. 

Cd( s ) has poles at s = 1,2, ...,q+ 1 in agreement with the general theory 
although in this case one cannot conclude that the asymptotic series for the heat 
kernel terminates. 

There are various pieces of information in the ( d (s) but we pass on to the 
vacuum energy. 

3. VACUUM ENERGIES 

The simplest cases are the even spheres. From (5) we see that ( d ( — 1/2) is zero 
because ( R (—2n) vanishes for n = 1,2, — Hence (2) gives 

E 2q = for S 2q . 

For odd spheres, S 2q+1 , substitution of (4) into (2) gives the finite value, 




in terms of Bernoulli numbers. 

As a function of dimension E d oscillates about zero. Numerical evaluation 
gives, for example, 

E 3 = l/240a 

E 5 = -31/60480a 

E 7 = 289/3628800a 

and a graphical plot can be found in fig.l, together with a hand interpolation. 

The simplicity of this calculation is due, in part, to the assumption of conformal 
coupling in d+ 1 dimensional space-time. For any choice of £ other than (1) it is a 
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little harder to determine the analytical structure of the zeta functions. A different 
£ is equivalent to giving the field a mass k and leaving £ as in (1). The eigenvalues 
(3) are then replaced by k 2 + (n + u — 1/2) 2 a~ 2 . All the algebraic manipulations 
go through exactly as before. The only difference in the final answer for the zeta 
functions Cd( s ) ^ n (4) an d (5) is that the Riemann zeta functions are replaced by 
zeta functions related to those investigated by Ghika and Visinescu (1978), Gibbons 
(1977) and Ford (1980) for example. 

For odd spheres a renormalisation is now necessary, but not for even spheres. 
(This is known from general considerations.) A more detailed discussion will be 
presented at another time. For example on S 2 for small Ka, E ~ — it 2 (kci) 4 / 256a 
while for large na, E ~ — («a) 2 /24a 

Another system for which there is no infinite renormalisation is the cuboid with 
Dirichlet or periodic boundary conditions. A summary of some results occurs in 
the next section. 

4. CUBIC DOMAINS. 

For periodic boundary conditions in a <i-cube of side a, the zeta function is 
given by 



in terms of the Epstein zeta function Z d . (We use the notation of Erdelyi 1953.) 
For d = 1, 2, 4, 6 and 8 'closed' expressions for Z d are known (e.g. Siegal 1944.) 



Z 1 (2s,l) = 2( R (2s) 
Z 2 (s,1) = Kr(s)P(s) 
Z 4 (s/2,l)=8(l-2 2 - 2s )( R (s)( R (s-l) 
Z 6 (s/3, 1) = 16( R (s - 2)p(s) - 4p(s - 2)( R (s) 
Z 8 (a/4, 1) = 16(1 - 2 l ~ s + 2 4 - 2s )( R (s) ( R (s - 3) 



permitting an easy evaluation of the vacuum energies, E d . We find the values for 
the pair (d,-aE%), (1,0.524), (2,0.719), (3,0.838), (4,0.932), (5,1.02), (6,1.122), 
(7, 1.23), (8, 1.40). We have included the three dimensional case computed using 
the lattice sums of Lennard-Jones and Ingham (1925). (See e.g. Unwin 1982). The 
d = 5 and 7 values were obtained by interpolation in lieu of a proper numerical 
evaluation. The graph is depicted in fig. 2. 
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For Dirichlet boundary conditions the zeta function is more complicated be- 
cause of the restriction on the modes but elementary manipulation of the summa- 
tions yields the general relation between cubic Dirichlet and periodic zeta functions, 

r=0 ^ ' 

and inversely 




From (2) the same relations hold for the vacuum energies. For example for the 
4-cube 

= 0.0063/a. 

The results are gathered on fig. 3 and again we see an oscillatory behaviour 
about zero as the dimension increases. 

Periodic boundary conditions turn the <i-cube into a d-torus, which, being 
multiply connected, can also support twisted fields, as is well known. 

So far as the modes go, the twisting amounts to replacing integers n by n — 1/2 
so that the fields reverse sign on traversing the appropriate circumference. In this 
case the torus zeta function is 

CJ«=(£) 2 *^o)(^,i) 

where the <i-dimensional vector a. has components equal to either 0, for an untwisted 
mode, or —1/2 for a twisted one. 

Cj(s) is to be evaluated at s = —1/2 and so the functional relation for the Z d 
must be used to make the argument positive to facilitate numerical computation. 
One finds 

which is the image form, as explained in Dowker and Banach (1978). 

The zeta function here has a denominator which is a sum of squares of integers, 
m 2 say, but the signs of the terms in the defining summation can be either +1 or 
— 1 since the numerator, exp(27rzo;.m), is just a product of —Is. If there are b 0s 
and a ( — l/2)s in a, (a + b = d), then Z d (0, a) is what Zucker (1974) calls S(a, b) 
and his expressions can be used to give the required numbers with no real work. 
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For simplicity we consider the case when all directions are twisted a = d, b = 0. 
Then 



Ej = -— r( a )5(d,0) 



,_ d+1 



with, corresponding to (5), 

S(l,0) = -2 V (s) = -2(l-2 1 - s )( R (s) 

S(2,0) = -± v (s)P(s) 

S(4,0) = -8 V (s-l)rj(s) 

5(6, 0) = -16 V (s - 2)(5{s) + A(5{s - 2)( R (s) 

5(8,0) = -16»7( a -3)C fl (a). 
Writing the results in the form (d,aE%) we find (1,0.262), (2,0.21), (3,0.196), 
(4,0.202), (6,0.266), (8,0.396). 

Again the d = 3 value has been evaluated using the numbers of Lennard-Jones 
and Ingham (1925) (5(3, 0)(s) = A'^ — A^ in their notation and = d(2s) in Zucker's 
(1975)). 

The graph (fig. 2) of these values shows a minimum at d = 3. Thus if the 
universe were toroidal and if the vacuum energy of twisted scalar fields dominated 
everything else, the universe would be prevented from sprouting extra toroidal di- 
mensions on energy grounds. 

This conclusion is too artificial to be significant but it is quite possible that, in a 
more realistic theory, if the dimension of the universe were a dynamical variable, and 
initially undecided, then vacuum energy considerations could force it preferentially 
to be four (or to be whatever the preferred dimension is.) 

It is straightforward to generalise Maxwell theory to d dimensions, and for 
perfect conductor boundary conditions in the <i-cube it is easy to show that the 
vacuum energy, Ef M is related to the scalar Dirichlet value by 

ES M = (d-l)ES + dESL 1 . 

The three dimensional case is discussed by Lukosz (1973a,b) and the general 
case by Ambj0rn and Wolfram (Cal.Tech. report CALT 68-855). (Incidentally the 
criticism in the footnote on p. 14 of this last reference is unfounded.) 

Substitution of the numbers gives values for (d, Ef M ), (2, -0.221), (3, 0.0915), 
(4,-0.045), (5,0.0211), (6,-0.01), (7,0.005). In view of the interpolation and of 
the occurrence of cancellations in formulae like (6) these numbers should not be 
taken to be very accurate. However we do see again the oscillatory behaviour with 
dimension. 
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5. DISCUSSION. 



The above detailed calculation can be extended in several ways. Firstly it 
would be interesting to consider deformed spheres. In the case of an S 3 with one 
scale differing from the other two, the neutrino and photon vacuum energies can 
be found exactly, up to a numerical integration. The same should be true for S d . 
Deformed cubes (cuboids) have already been discussed (Dowker and Banach 1978, 
Mamaev and Trunov 1979, Ambj0rn and Wolfram loc.cit., Unwin 1982). 

Secondly, higher spins can be considered and it would be interesting in this 
case to see whether the heat kernel asymptotic series terminates. In fact it is an 
interesting problem to determine the exact conditions for such a termination in an 
arbitrary curved space and to elucidate the relation between the termination and 
the WKB approximation. 

Thirdly, everything can be done at finite temperature. This is simply a tech- 
nical problem and easily amenable to calculation. It is relevant to bag and early 
universe theory. 
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FIGURE CAPTIONS 



FIG 1. 

Vacuum energy of a conformally coupled scalar field on the <i-sphere plotted 
against d. The non-integer interpolation is purely cosmetic. 
FIG 2. Vacuum energies on the <i-torus for twisted and untwisted scalar fields. 
FIG 3. Vacuum energy in the <i-cube for Dirichlet boundary conditions. 
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This figure "figsl, 2.jpg" is available in "jpg" format from: 



http://arxiv.org/ps/! 106. 3657vl 



This figure "fig3.jpg" is available in "jpg" format from: 



http://arxiv.org/ps/! 106. 3657vl 



